Diffractive resonant radiation emitted by spatial solitons in waveguide arrays: towards 

a spatiotemporal supercontinuum generation 
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We study analytically and numerically a new kind of diffractive resonant radiation emitted by 
spatial solitons, which is generated in waveguide arrays with Kerr nonlinearity. The phase matching 
condition between soliton and radiation is derived and agrees well with direct pulse propagation 
simulations. The folded dispersion due to the Brillouin zone leads to a peculiar anomalous soliton 
recoil that we describe in detail. A linear potential applied across the array generates the analogue 
of the Raman self- frequency shift in optical fibers, only now applied to the wavenumber. We 
demonstrate that it is possible to mimic closely temporal fiber-optical dynamics, unveiling the 
new effects of wavenumber- supercontinuum generation and the compensation of the 'soliton self- 
wavenumber shift' by the emitted diffractive radiation. This work paves the way for designing unique 
optical devices that generate spectrally broad supercontinua with a controllable directionality. 



Introduction — Waveguide arrays (WAs) consisting 
of identical, equally spaced waveguides, present a unique 
discrete platform to explore many interesting fundamen- 
tal phenomena such as discrete diffraction [ij, discrete 
solitons 0, m, and photonic Bloch oscillations [4-7^. 
In applications, 2D networks of nonlinear waveguides 
with discrete solitons may be useful for designing signal- 
processing circuits [8]. Recently, WAs have also been 
used to mimic relativistic phenomena typical of quantum 
field theory, such as Zitterbewegung [9], Klein paradox 
K 



fermion pair production [11] and the Dirac equation 



The concept of dispersive resonant radiation (DisRR), 
which emerges due to higher-order dispersion (HOD) 
terms, has been well studied in the last decade in the 
temporal case for optical fibers jl5l4l9( . When an ul- 
trashort pulse is launched into optical fibers, a DisRR 
due to the phase matching between the fiber and the 
soliton group velocity dispersion (GVD) generates one 
or more new frequencies [3, [l^ . This DisRR, together 
with other nonlinear effects such as self- and cross-phase 
modulation, soliton fission [13], and stimulated Raman 
scattering [14], are the main ingredients of the so-called 
supercontinuum generation (SCG) [lO, HH, particularly 
in highly nonlinear photonic crystal fibers [22]. SCG is 
among the most important phenomena in nonlinear fiber 
optics which has led to a number of important technolog- 
ical advances in various fields, such as spectroscopy and 
medical imaging [23], metrology [24], and the realization 
of broadband sources [25]. 

Inspired by advances in DisRR studies in the last 
decade, in this Letter we prove theoretically the existence 
of a new kind of resonant radiation - namely the diffrac- 
tive resonant radiation (DifRR), which occurs when a 
continuous- wave (CW) beam or a relatively long pulse is 
launched into WAs. Similarities and differences between 
DisRRs and DifRRs are analyzed. We show that when 
the phase matching condition is satisfied, a spatial soli- 
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FIG. 1: (Color online) (a) Solid blue line: WA dispersion Kz 
vs. K,. Dashed red line: D2 vs. k,, showing the two zero- 
diffraction points located at At = ±7r/2. (b) Wavenumber 
krr of the generated DifRR, as a function of the input soli- 
ton wavenumber kq. In both (a,b) the gray shaded area in- 
dicates the region where bright solitons can propagate. (c,d) 
Beam propagation in the (n, ^)-plane (c) and {k, z)-plane (d). 
Parameters are Aq = 0.8, C = 1.2, kq = 0.7, N = 100. 



ton emits DifRR with a new well-defined direction, i.e. 
transverse wavenumber. Moreover, due to the periodic- 
ity of discrete systems, and thus the existence of a Bril- 
louin zone, many unusual effects which cannot exist in 
continuous media can now occur, such as the anomalous 
solitonic recoil and the anomalous compensation of soli- 
ton self- wavenumber shift, described in this Letter for the 
first time. This paves the way for producing controlled 
spatiotemporal supercontinua in both the frequency and 
wavenumber domain, i.e. broadband radiation with a 
controlled directionality. 

Phase matching condition for the diffractive resonant 
radiation — Light propagation in a discrete, periodic 
array of Kerr nonlinear waveguides can be described, in 
the CW regime, by the following dimensionless equation 
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arrives at the following equation: 



^ ^^"dz^ +c[an+i(2;) + an-i(2;)] + \an{z)\^ an{z) = 0, (1) 

where is the electric field amplitude in the nth 
waveguide, z is the longitudinal spatial coordinate, c = 
C/{'yP)^ where C is the coupling coefficient (in units of 
a wavenumber) resulting from the field overlap between 
neighboring waveguides, 7 is the nonlinear coefficient of 
a single waveguide, and P is the input beam peak power. 
By using the stationary discrete plane wave solution for 
the n-th waveguide an{z) = ao ex.-p[i{nkxd -\- nzz)] one 
arrives, in the linear case, at the dispersion relation be- 
tween and kx 



i^zikx) = 2ccos{kxd), 



(2) 



where d is the center-to-center spacing between two ad- 
jacent waveguides, and kx is the transverse wavenum- 
ber, see solid blue line of Fig. [IJa). It is clear from 
Eq. ([2j) that is periodic in a^: = kxd, which represents 
the phase difference between adjacent waveguides. Thus, 
within the coupled mode approximation, it suffices to in- 
vestigate the first Brillouin zone of the folded dispersion, 

— TT < K < TT. 

Since a typical input beam has a finite width cover- 
ing several waveguides, its Fourier spectrum has a cer- 
tain bandwidth with a central transverse wavenumber 
/^o, which is fixed by the input angle of incidence of the 
exciting beam. We can then use a Taylor expansion of 
Eq. ^ as follows: 
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where Ah 



and Dj^ 



the mth order diffr active Taylor coefficient [thus, Di = 
—2csm{f<io) ^ D2 = — 2ccos(ai:o), etc.]. In Fig. [TJa) we 
plot a typical curve for 1^2 (^) (dashed red line), show- 
ing the existence of two zero-diffraction points located at 
hz = ±7r/2. In a sense, this shape of D2 is analogous to 
the GVD of photonic crystal fibers in the temporal case 

122.]. 

Following Refs. [71, l26|, we now approximate the dis- 
crete variable n with a continuous one. This is jus- 
tified since we shall use pulses and solitons that ex- 
tend for several waveguides. Defining now n as a 
continuous variable of the distributed amplitude func- 
tion ^(n,z) = an,2 exp(— z/^o^), we eliminate the zero- 
th order term /^^(^o) which is responsible for a gen- 
eral phase evolution through the substitution ^(n,z) 
"^{n, z) ex.p[ihZz{i^o)z]. The first order term, —iDidn, 
takes into account the transverse velocity and can also 
be eliminated by introducing a comoving frame, n 
n -\- Diz. After dropping these two low-order terms one 
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*(n,2) = 0. 
(4) 



Equation (|4]) is formally identical to the well-known 
generalized nonlinear Schrodinger equation (GNLSE), 
which describes the evolution of pulses in a single op- 
tical fiber, plus HOD terms [25^. In Eq. (j4]) we have 
the transverse spatial variable n instead of the temporal 
variable t of the conventional GNLSE. Unlike the tempo- 
ral GNLSE, where a Taylor series for the fiber dispersion 
can usually be expanded up to a small number of terms 
(because HOD coefficients become rapidly very small), 
in Eqs. (|3]) and (|4]) many higher-order diffraction terms 
Dm>2 should be taken into account, since their absolute 
values will be either |2csin(Ai:o)|, or |2ccos(/>:o)|, and the 
sum only converges due to the factorial in the denomina- 
tor. 

In the temporal version of the GNLSE, it is well-known 
that a temporal soliton propagating in a fiber emits 
small-amplitude, dispersive and quasi-monochromatic 
waves at well-defined frequencies (the DisRR) when the 
linear fiber dispersion and the nonlinear soliton disper- 
sion are matched [Til, [HI. It is thus natural to con- 
jecture that, in a WA, a spatial soliton, which in the 
continuous variables approximation extends over several 
waveguides, emits during the propagation a similar kind 
of small-amplitude diffractive radiation, within a narrow 
wavenumber range, due to the phase-matching between 
the spatial soliton nonlinear dispersion and the linear 
array dispersion given by Eq. (|2]). By using the per- 
turbation approach which was developed for DisRRs in 
[18], here we derive the phase-matching condition for the 
DifRR in a similar way. We first find the unperturbed 
soliton solution of Eq. (|4]) where all diffractive terms 
Dm>3 are dropped. The soliton solution of Eq. (j4]) is: 



aso\{z^n) = Aosech 
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exp{ikso\z), (5) 



where /cgoi = ^o/^ spatial soliton longitudinal 

wavenumber, identical to its temporal counterpart. The 
bright soliton solution ^ only exists when 2ccos(/t:o) > 0, 
i.e. only in half of the Brillouin zone, where — 7r/2 < 
^0 < 7r/2. Now we look for the linearized dispersion re- 
lation of plane wave solutions of Eq. (|4]) , by substituting 
ex.-p[i{k\inZ + AtvTi)] into Eq. (|4]) and using Eq. (|3j). We 
obtain: 

k\in(AK,) = — ^Atv^ = 2c\cos{K,)—cos(K,o)-\-s'm{K,o)AK] 
^-^ ml 

m>2 

(6) 

In Eq. (|6|), /^o is the central wavenumber (which is related 
to the incident angle) of the incident beam, while Ak is 
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the detuning from /^o, and k = /^o + A/^. Energy exchange 
between radiation and sohton is possible for those values 
of Ahi that satisfy 



(7) 



where ksoi is constant and has been defined above. This 
phase matching condition, an implicit equation for the ra- 
diation wavenumber detuning A/^, is the central result of 
this Letter, and accurately predicts the DifRR wavenum- 
ber as we shall show soon. If we follow what is commonly 
done for optical fibers, and we take into account only D2 
and Ds in Eq. ([6]), and ignoring the power dependence 
(^soi 0)5 one can easily get the approximate DifRR 
wavenumber in the form /^rr k,'^^ = a^o + 3/tan(Ai:o). 
Such approximations are perfectly fine in fiber optics 
when dealing with DisRR - they lead to very accurate 
predictions of the DisRR frequency. However, the same 
approximation is not good enough for the case of WAs, 
since, as explained above, the coefficients {Dm/Tn\) de- 
cay not as rapidly as in the temporal case, and a large 
number of orders must be taken into account. However, 
even if not explicit as in the case of the temporal DisRR, 
Eq. ([7]) is exact and can be easily solved numerically. 
Emission of DifRR and soliton anomalous recoil — 
We now prove numerically the formation of DifRR 
and the accuracy of the predictions made by the phase- 
matching condition Eq. ([7]). In Fig. [TJb) we show the 
DifRR wavenumber /^rr = /^q + A/^ as a function of the 
input soliton wavenumber (which is related to the an- 
gle of incidence) hcq. The blue solid curve is obtained 
by finding numerically the roots A Hi of Eq. ([7]), while 
the dashed red curve shows the approximated analyti- 
cal expression given in the previous section. It is clear 
that K^'^iYi is not accurate enough to be used in practice, 
when compared to the solid line. In Fig. [Hb) we depict 
the full range of the first Brillouin zone for completeness, 
but only the interval — 7r/2 < < 7r/2 (indicated by 
a gray shaded area), in which pulses experience anoma- 
lous diffraction, should be considered. Parameters used 
in Fig. [TJare: ^40=0.8; c=1.2; /^o=0-7. For these param- 
eters, in the range 0.235 < \i<io\ < 7r/2, one can find only 
one solution for /^rr, but when < \hZo\ ^ 0.235, Eq. 
([7j) shows several roots [see the solid blue curve in Fig. 
[T](b)]. Thus, one should expect to simultaneously gen- 
erate several DifRRs with different wavenumbers in the 
latter interval. However, full numerical simulations of 
Eq. ^ show that only the solution corresponding to the 
branch that is the closest to the central horizontal axis 
(i.e., the axis hz^iR = 0) can be generated and observed, 
and all other DifRR waves corresponding to roots from 
other branches are too weak to be seen numerically, since 
the overlap between the soliton spectral tail and the radi- 
ation wavenumbers becomes exponentially small. When 
hzo = 0, i.e. for a direct incidence of the input CW beam, 
there is no solution for Eq. (O, regardless of the parame- 
ters used. This is also confirmed by the direct simulation 



ofEq. ©. 

The evolution of a CW beam along z according to 
Eq. ([1]) is shown in Fig. ttl^c), for an input beam 
<^in(^) = Aosech[nAo / ^/2ccos{Hio)]e'^'^°'^ (i.e. the approx- 
imate soliton solution in the continuous limit), and for a 
WA made of = 100 waveguides. After some propaga- 
tion, around z 3, a DifRR is emitted by the soliton. 
The evolution of the Fourier transform of the field a{n) 
of Fig. [TJc) along z is shown in Fig. [IJd). The dashed 
green horizontal line represents the input wavenumber 
{hzo = 0.7), while the solid white line is obtained by 
solving Eq. ([7]) numerically, showing excellent agreement 
with the pulse propagation. In Fig. [Hb) one can notice 
that the soliton emits the DifRR with a positive detuning 
Ahi when < hcq < 7r/2. For instance, when a^o = 0.7, 
then from Eq. ([7]) one gets tZRR = /^o + A/^ ~ 3.53. 
However, since the Brillouin zone has a limited exten- 
sion, when 27T > krr > tt the DifRR will be emitted 
with a negative detuning due to the folding of the band 
structure. In the example shown in Figdjd) the effective 
DifRR wavenumber will be equal to hZRR — 27r o:^ —2.75 
[see the white solid line in Figdjd)]. This means that 
in real space the soliton, instead of recoiling in an op- 
posite direction than the DifRR, will recoil towards the 
DifRR itself, see the white arrows in Fig. [U^c). The same 
phenomenon occurs in the wavenumber space: the soli- 
ton spectral momentum, instead of recoiling away from 
the radiation, moves slightly towards it [see Fig. [Hd)]. 
We call this unique effect anomalous recoil. The final 
spectrum of Fig. [Hd) can thus be best described as su- 
percontinuum generation in wavenumber space. 
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FIG. 2: (Color online) (a,b) PBOs in the linear regime, shown 
in both the (n, 2;)-plane (a) and the {k, z)-plane (b). Param- 
eters are c = 1, a^o = 0, = 0.8, a — 0.5 and N — 35. (c,d) 
Anomalous compensation of SSWS in the nonlinear regime 
depicted in the (n, ^)-plane (c) and {k, z) -plane (d). Param- 
eters are c = 1, /t:o = 0, = 0.6, a — 0.01 and N — 300. 
White solid lines in (b,d) represent the wavenumber calcu- 
lated by using the formula Hi(z) = kq -\- az, making sure that 
K G [— TT, tt]. 
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Compensation of the soliton self-wavenumber 
The Raman soliton self-frequency shift (SSFS) 
is a fundamental effect acting on temporally shortj 
in silica fibers, and plays a crucial role in SCG [2| 



)ulses 
1- In 

fibers, SSFS leads to a continuous shift of the central 
frequency ujq of the soliton towards the red part of the 
spectrum [27]. When ujq is sufficiently close to one of the 
zero-GVD point of the fiber, DisRR is emitted, and the 
SSFS can be compensated [li,[29|. 

Here we show that an analogue of the Raman SSFS 
can be obtained in WAs by applying an external lin- 
ear potential across the transverse coordinate n, which 
changes the propagation constant along the array in a 
linear fashion, by using for instance the electro-optic [3j 
or thermal-optic effect (Hf. In the continuous limit, this 
results in adding a term anan{z) to the the left-hand 
side of Eq. ([T]), and thus a new term an"^ to the left- 
hand side of Eq. (j4]), with a depending on the spe- 
cific process used to create the potential. By using the 
same moment method developed for calculating the shift 
rate of SSFS in the temporal case 25, |30|, the soliton 
self-wavenumber shift (SSWS) rate turns out to satisfy 
the simple expressions for the central wavenumber and 
the central position of the soliton, hz{z) = tzo -\- az and 
no{z) = no{0) -\- ccos{hzo){az'^ -\-2hZoz) ^ thus the wavenum- 
ber shift rate being simply equal to a, which is also pro- 
portional to the rate of acceleration or deceleration [for 
sgn(a) = ±1 respectively] in the n space. 

Actually the above formula of hz{z) is exact in the ab- 
sence of Kerr nonlinearity. One immediate effect of the 
wavenumber shift in the linear regime is the onset of the 
well-known photonic Bloch oscillations (PBOs), see Refs. 
[J^]. We can interpret PBOs in the following way. In 
Fig. Efb) a linear beam is launched in a WA made of 
iV = 35 waveguides at normal incidence (/^o = 0), but due 
to the wavenumber shift described above, its wavenum- 
ber changes linearly towards positive detunings Ak, (for 
a > 0), and therefore it accelerates towards positive val- 
ues of no, as shown in Fig. [2fa). However, due to the fold- 
ing of the band structure, when the wavenumber reaches 
the value /t: = tt, it folds back into the first Brillouin zone, 
and thus becomes negative, see white solid lines in Fig. 
[2jb). This means that the beam must change direction, 
and starts decelerating towards negative values of no , and 
this process is repeated infinitely in a sinusoidal fashion 
with a period zq = 27r/a, see Fig. EJa). Therefore, the 
above formula for no(^) is only valid at the beginning of 
propagation, before the occurrence of the folding effect. 

The analogy between SSFS and SSWS appears instead 
in the nonlinear regime. Even though the soliton input 
wavenumber is far from the zero diffraction points ([lo- 
cated at ±7r/2 in the dispersion of Fig. [U^a)], eventually 
the SSWS will push the soliton wavenumber in proximity 
of one of these two points, depending on the sign of a. 
At this moment, one can conjecture that the soliton will 
start emitting a strong DifRR that is able to compen- 



sate the SSWS. This is indeed the case, as shown in Figs. 
[2](b,c). In Fig. [2fc) we show the z evolution of a soliton 
in the n space with an initial position no(0) = 0, in a 
potential with a = 0.5. We use a WA made of = 300 
waveguides, but this is done only in order to make all the 
figures as clear as possible for the reader, while the con- 
cept works for an arbitrary number of waveguides. The 
soliton starts accelerating (for a > 0) towards positive 
values of no or decelerating (for a < 0) towards nega- 
tive values of no continuously, until at a specific moment 
(around z 70) it emits the DifRR. In Fig. Efd) we show 
the same as Fig. [2fc), but in the wavenumber space. The 
soliton starts with zero central wavenumber (i.e. normal 
incidence, hzo = 0), shifts continuously and linearly its 
wavenumber (the white solid line shows the predicted 
shift, hi = hiQ -\- az^ for the central wavenumber) until a 
considerable amount of soliton energy reaches the zero 
diffraction point = 7r/2, which triggers the emission 
of a strong DifRR. This emission is so strong that the 
wavenumber shift of the soliton is compensated, a pro- 
cess that is an exact analogue of the SSFS compensation 
in fibers pj|. 

However, note that in Fig. [2fc,d), an anomalous re- 
coil identical to that of Figs. [Hc,d) occurs: the DifRR 
should be emitted by the soliton in the positive detunings 
Ak, > 0^ but the folding of the Brillouin zone makes the 
radiation appearing in the negative detunings - and thus 
in this case also the soliton central wavenumber shifts to- 
wards the DifRR, and not away from it, as it happens in 
optical fibers. 

Conclusions — We have discovered the existence and 
studied the properties of diffractive resonant radiation 
emitted by discrete spatial solitons in waveguide arrays. 
Due to the periodicity, several effects such as an anoma- 
lous soliton recoil and the anomalous compensation of the 
soliton self-wavenumber shift, which have no counterpart 
in continuous systems, are revealed. The analysis of these 
unusual phenomena could be useful for understanding su- 
percontinuum generation and many other spatiotemporal 
effects in waveguide arrays, and are applicable to virtu- 
ally any nonlinear discrete periodic system supporting 
solitons, therefore making our results very general and of 
relevance for a number of very diverse communities. 

We acknowledge Ulf Peschel for useful discussions. 
This work is supported by the German Max Planck So- 
ciety for the Advancement of Science (MPG). 
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